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Understanding nuclear effects in parton distribution functions (PDF) is an essential component 
needed to determine the strange and anti-strange quark contributions in the proton. In addition 
Nuclear Parton Distribution Functions (NPDF) are critically important for any collider experi- 
ment with nuclei (e.g. RHIC, ALICE). Here two next-to-leading order ^-analyses of NPDF are 
presented. The first uses neutral current charged-lepton [l^A) Deeply Inelastic Scattering (DIS) 
and Drell-Yan data for several nuclear targets and the second uses neutrino-nucleon DIS data. We 
compare the nuclear corrections factors (F^'^/Fj^) for the charged-lepton data with other results 
from the literature. In particular, we compare and contrast fits based upon the charged-lepton DIS 
data with those using neutrino-nucleon DIS data. 
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1. Introduction 

Parton distribution functions (PDFs) are extremely important in high energy physics as they 
are needed for the computation of reactions involving hadrons. For this reason various groups 
present and update precise global analyses of PDFs for protons [[T], |^ |3|, ^ and nuclei [||, |^ 
PDFs are non-perturbative objects determined by experimental input and the data come from dif- 
ferent processes such as the Deep Inelastic Scattering (DIS), Drell-Yan (DY) and jet production. A 
lot of data is extracted from free protons but there is a large share which comes from analysis of 
data on nuclear targets. Most prominently, the neutrino DIS on heavy nuclei is proving to be very 
important for precise determination of the flavor components of the PDFs and it gives the most 
precise information on the strange quark PDF. The knowledge of the strange quark PDF has an 
impact on the precision of W and Z boson measurements at the LHC. Data taken on nuclear targets 
are included in the proton analysis using the nuclear correction factors which are very often based 
on a specific model [^]. The other option is to use the data on nuclear targets and to extract from 
them the nuclear parton distribution functions (NPDFs) in order to construct the nuclear correction 
factors based on experimental input. NPDFs are also used in predictions for colUsions of nuclei at 
RHIC or at the LHC. 

Here, we present a new framework for a global analysis of nuclear parton distribution functions at 
next-to-leading-order (NLO). Then we use it to analyze the apparent differences between nuclear 
correction factors (Fj^^'/Fj^^'^) coming from charged lepton data and from neutrino DIS data. 

2. NPDF global analysis framework 

We introduce the global analysis framework and the analysis of NPDFs using the charged 
lepton DIS and Drell-Yan data for a variety of nuclear targets. The analysis is performed based on 
the same principle as the proton analysis of ||^. The input distributions are parameterized as 

xfk{x,Qo) = C(ix'''{\-xy^e''^''{\+e''''xY' k = u^,dy,g,u + d,s,s, (2.1) 

dix,Qo)/uix,Qo) = cox^'(l-#' + (l+C3^)(l-#^ 

at the scale Qo= 1.3 GeV. The different nuclear target materials are treated by introducing a nuclear 
A-dependence in the Ck coefficients: 

Ck^Ck{A) = Ck,o + Ck,i{l-A-'''''), k = {l,...,5}. (2.2) 

The advantage of this construction is that in the limit A — > 1 we recover the original proton parame- 
terization with Ck,o as the coefficients. Using this framework, we construct a global fit to the charged 
lepton DIS data and Drell-Yan data (for details see [0]). The coefficients ctfl were based on the 
results of the proton global fit presented in [ [To| ] where the influence of nuclear targets on proton 
PDFs was minimal. In the analysis, we have applied standard kinematic cuts of 2^ „, =2.0 GeV, and 
Wcut = 3.5 GeV. Performing the global fit to the data (708 data points after the cuts are applied), 
we obtain an overall ;^^/dof of 0.946 with about 32 free parameters. The results of the global 
fit, the coefficients Ck,i and cic,2, give the A-dependence of the generalized coefficients Ck{A) (see 
Fig. |l|) and these coefficients determine the parton distribution functions for bound partons inside 
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Figure 1: The A -dependent coefficients cii{A), k — {1,5}, for the up-valence (a) and down-valence PDF (b) 
as a function of the nuclear A. The different coefficients Cjt(A) correspond to following lines: ci- solid (red) 
line, C2- long dashed (blue) line, C3- dashed (green) line, C4- dash-dotted (magenta) line, C5- dotted (brown) 
line. 



a nucleus (see Fig. The nuclear effects are typically given in terms of nuclear correction factors 
R[F2^^''] = /Fj'^'^ and we show the nuclear correction factors resulting from our fit to charge 
lepton data in Fig. |3[ The nuclear structure functions and ^2^^'" are both defined according to 

^■^^•^^ {x. Q) = \ Ffl^ {x, Q) + F^" (x, 2) . (2.3) 

These structure functions can be computed at next-to-leading order as convolutions of the (nuclear) 
PDFs with the conventional Wilson coefficients, i.e., generically 

F;^^^\x,Q) =l^Q,® fl^^^K (2.4) 



The difference between ^2^^ and f!^^''^ is in using different PDFs {F2'' uses nuclear PDFs and F2 



Fefi 



uses proton PDFs) in Eq. ( |2.4| ) 



3. Nuclear correction factors from neutrino DIS 

The result of an analysis of NuTeV neutrino DIS cross-section data performed in [ p^ showed 
a deviation from the standard result of the analysis of charged lepton DIS and DY data. This 
can be clearly seen when comparing the different nuclear correction factors in Fig. ^. As the 
different nuclear correction factors were not obtained in completely identical frameworks, we first 
used the NPDF framework introduced in the previous section and in JTl] ] to re-analyze the NuTeV 
neutrino DIS data. Using the same kinematic cuts, we obtain a fit to 2310 neutrino and anti-neutrino 
cross-section and di-muon data points. The result of the fit (see Fig. ^) is in agreement with the 
previous analysis and confirms the difference between the nuclear correction factors mainly in the 
intermediate x-region. The obvious difference poses a question if a compromise nuclear correction 
factor can be found which would accommodate both the charged lepton and neutrino data. In order 
to construct a compromise fit we use all data used for the charged lepton fit (708 data points), 
NuTeV and Chorus neutrino and anti-neutrino cross-section data and NuTeV and CCFR di-muon 
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Figure 2: We display the (a) xmv(x), (b) xdv{x), (c) xg{x) and (d) xs{x), PDFs for a selection of nuclear A 
values ranging from A = {1,207}. We choose Qq = 1.3GeV. The different curves depict the PDFs of nuclei 
with the following atomic numbers (from top to bottom in (c) atx = 0.01) A = 1,2,4,8,20,54, and 207. 



data (3134 data points). To avoid the neutrino data to dominate the analysis, we apply a weight 
factor 1/2 to the coming from the neutrino and anti-neutrino cross-section data. The result of 
such a compromise fit with the weight 1/2 is shown in Fig. Although the nuclear correction 
factor from the compromise fit seems to be compatible with the charged lepton data, to draw a 
firm conclusion a further investigation is necessary. A detailed analysis is postponed to the next 
publication (see also analysis in |[T3[]). 



4. Conclusions 

We presented a framework for a global analysis of NPDFs at next-to-leading order QCD 
closely linked to a proton analysis. We used this framework to analyze the discrepancies between 
the nuclear correction factors stemming from the analysis of charged lepton DIS and DY data and 
the ones coming from neutrino DIS data. We confirm the differences found in a previous analysis 
and we presented preliminary results on compromise fit combining charged lepton and neutrino 
data. A much more detailed analysis is postponed to a later publication. 
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Figure 3: The nuclear correction factor as a function of .n: for — 20GeV^. Figure-(a) shows the fit (fit 
C) using charged-lepton-nucleus DIS and DY data whereas Figure-(b) shows the fit using neutrino-nucleus 
data (fit A2 from Ref. [|l2|]). Both fits are compared with the SLAC/NMC parameterization, as well as fits 
from Kulagin-Petti (KP) (Ref. g]) and Hirai et al. (HKN07), (Ref. [g). 
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Figure 4: As above but here Figure-(a) shows a compromise fit using all charged lepton and neutrino data 
with neutrino cross-section data weighted down by a factor of 1 /2 whereas Figure-(b) shows the fit using 
only neutrino and anti-neutrino DIS data. 
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